Coulomb interaction in graphene: Relaxation rates and transport 
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We analyze the inelastic electron-electron scattering in undoped graphene within the Keldysh 
diagrammatic approach. We demonstrate that finite temperature strongly affects the screening 
properties of graphene, which, in turn, influences the inelastic scattering rates as compared to 
the zero-temperature case. Focussing on the clean regime, we calculate the quantum scattering 
rate which is relevant for dephasing of interference processes. We identify an hierarchy of regimes 
arising due to the interplay of a plasmon enhancement of the scattering and finite-temperature 
screening of the interaction. We further address the energy relaxation and transport scattering 
rates in graphene. We find a non-monotonic energy dependence of the inelastic relaxation rates in 
clean graphene which is attributed to the resonant excitation of plasmons. Finally, we discuss the 
temperature dependence of the conductivity at the Dirac point in the presence of both interaction 
and disorder. Our results complement the kinetic-equation and hydrodynamic approaches for the 
collision-limited conductivity of clean graphene and can be generalized to the treatment of physics 
of inelastic processes in strongly non-equilibrium setups. 



I. INTRODUCTION 

Graphene!'- is a two-dimensional (2D) material with 
a quasi-relativistic dispersion law that has attracted an 
outstanding attention of leading experimental as well 
as theoretical groups all over the world. In 2004, re- 
searchers at Manchester University first succeeded in ex- 
perimental isolation of a monoatomic graphite layer — 
graphene — on an insulating substrate 3 -^. This technolog- 
ical breakthrough was immediately followed by transport 
measurements^* 6 - which have shown remarkable proper- 
ties related to Dirac nature of the charge carriers. In 
particular, a short and wide sample of clean graphene 
exhibits a pseudo-diffusive charge transport 7 , with the 
counting statistics equivalent to that of a diffusive wire 8 -! 9 -. 
This equivalence has been confirmed in recent measure- 
ments of conductance and noise in ballistic graphene 
flake a 10 ' 11 . In contrast to conventional metals, ballistic 
graphene near the Dirac point conducts better when po- 
tential impurities are addec— ~— . 

Quantum interference in disordered graphene is also 
highly peculiar due to Dirac nature of carriers. In par- 
ticular, at the Dirac point, the minimal conductivity^* 6 - 
~ e 2 /h is "protected" from quantum localization in the 
absence of intervalley scattering^ 6 - or in the case of a 
chiral-symmetric disorde r 17 ' 22 . Away from Dirac point, 
the concentration dependence of graphene conductiv- 
ity in diffusive samples depends strongly on the nature 
of scatterers— The experimentally observed (approxi- 
mately linear) dependence in most of the samples may 
be explained by strong impurities creating resonances 
near the Dirac point ("midgap states") 17 ' 18 , yielding 
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dominant type (or types) of disorder and the correspond- 
ing disorder strength depend, of course, on technology of 
the sample preparation. 

How important is the electron-electron interaction in 
graphene? Recent works demonstrated manifestation 
of the interaction in dephasing rate providing the cut- 
off to quantum interference phenomena—"—, as well as 
in renormalization of conductivit y 26 ' 27 . Recent experi- 
ments also showed that the interaction plays a particu- 
larly prominent role in suspended graphene samples ! 28 ' 29 
In such samples the splitting of integer quantum Hall 
transition (attributed to interaction-induced spin/ valley 
symmetry breaking) is observed at magnetic fields as low 
as 2T (i.e. an order of magnitude less than in graphene 
structures on a substrate). Furthermore, one observes 
also fractional quantum Hall plateaus, which indicates 
the importance of electron correlations—"— Continuous 
advances in fabrication of high-quality graphene samples 
are expected to lead to further enhancement of the role 
of interaction in graphene. 

On the theoretical side, interactions may have a dra- 
matic impact on quantum electronic transport, especially 
in systems of reduced dimensionality" 33 . Very generally, 
the interaction-induced phenomena may be subdivided in 
two big classes, related to effects of renormalization and 
inelastic scattering, respectively. In the case of graphene, 
the interaction physics becomes an even more complex 
problem, in view of the "relativistic" dispersion of carri- 
ers. Interaction phenomena are particularly strong near 
the Dirac point, where the density of states vanishes (for 
a clean system), the screening by intrinsic carriers be- 
comes very inefficient, and the Drude conductivity (in 
the presence of disorder) is of the order of the conduc- 
tance quantum. 
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The dimensionless bare coupling constant a g = e 2 /vf 
describing the Coulomb repulsion in graphene samples on 
a SiC-2 insulating substrate is estimated to be 0.6 — 0.8 
and can be yet larger, a g ~ 2.2, in suspended graphene 
sheets. Transport experiments on the interaction correc- 
tion to conductivity in graphene on a substrate^ 7 - have 
found much smaller values of the interaction constant 
(or, equivalently, the interaction parameter r s ) which 
can be attributed to the renormalization effects which 
we briefly overview below. Furthermore, the value of the 
effective fine structure constant of freestanding graphene 
as inferred from the results of recent experiments on X- 
ray scattering in graphite^ was also found to be much 
smaller than the bare value of a. 

The interaction effects in the clean graphene have 
been considered in Rcfs. I35l - l4ll within the weak-coupling 
renormalization-group (RG) scheme justified for a large 
number of flavors of Dirac fermions. The main result 
of this consideratio n 35 ^ 6 is the renormalization of the 
Fermi velocity and hence of the interaction parameter. 
This perturbative renormalization group for graphene 
thus shows that in the clean case the Coulomb interaction 
is marginally irrelevant. In the disordered unified 
ballistic RG emerges describing renormalization of dis- 
order couplings and of the interaction 3 - 7 - - — . In Ref. |4l| 
the corresponding one-loop RG equations are derived for 
time-reversal-invariant disorder and in the limit of large 
number of valleys. 

At sufficiently strong a g , the Coulomb interaction in 
a clean graphene has been argued to give rise to vari- 
ous instabilities^ 3 - - — , in particular, to opening the gap 
in the Dirac spectrum of the quasiparticles (spontaneous 
mass generation). Physically, this instability leads to a 
phase transition from the semimetallic state to an insula- 
tor state, known as excitonic semimetal-insulator transi- 
tion* 3 -. Recently, this type of instability has been studied 
by effective mean-field-type approac h 43 i 44 , renormaliza- 
tion group method^ 6 -, and lattice simulation o 47 ' 48 . The 
mean-field consideration of the excitonic-type instability 
in graphene predicts a certain value of the coupling con- 
stant a g ~ 2 (close to that in a suspended graphene) at 
which the instability occurs. The predicted semimetal- 
insulator transition, however, has not yet been observed 
in experiments in zero magnetic field. One of the possi- 
ble reasons for that might be the presence of disorder—. 
In strong magnetic field — in the quantum Hall effect 
regime — the repulsive interaction between electrons may 
result in the Stoner instability 5 ^— giving rise to spon- 
taneous breaking of spin and/or valley symmetry. Ex- 
periments do show splitting of quantum Hall plateaus in 
strong magnetic fields^ 3 -, that is attributed to interaction 
effects. 

It has been debated in the literature whether the 
Coulomb interactions in graphene can be theoretically 
addressed within the standard Fermi-liquid-type pertur- 
bation theory that requires, in particular that the particle 
energy is much larger than the decay rate (r <C e). In 
order to understand to what extent graphene is a Fermi 



liquid, one has to explore interaction-induced inelastic 
collision rates. The inelastic quantum scattering rate at 
zero temperature (and finite quasiparticle energy) has 
been considered in Refs. I54ll55l It was found in these 
works that the behavior at Dirac point is rather peculiar 
and requires a careful incorporation of screening. 

Another highly nontrivial feature of graphene is that 
inelastic electron-electron collisions may limit the con- 
ductivity at the Dirac point without any disorder or 
phonon scatterin g 5 ^ 57 ' 59 . This peculiarity of graphene — 
which should be contrasted to conventional systems 
where interactions do not lead (in the absence of Umk- 
lapp scattering) to finite resistivity — is a consequence of 
the particle-hole symmetry and decoupling between ve- 
locity and momentum. As a result, although the total 
momentum of interacting particles is conserved during 
inelastic collisions, the total current may relax. The 
collision-limited conductivity of undoped graphene is 
found to be inversely proportional to a 2 and depends 
on temperature only through the renormalization of a g . 
It was found that the energy relaxation caused by in- 
elastic processes in graphene is fast, which allows one to 
treat the problem by using a relativistic hydrodynamic 
approac h 58 ' 60 ^ . 

In this paper we analyze the inelastic electron-electron 
scattering in graphene within the Keldysh diagrammatic 
approach. While we focus on the equilibrium situation 
in this work, we have in mind to extend the treatment of 
physics of inelastic processes to strongly non-equilibrium 
setups, which explains why we prefer to work in the 
framework of the Keldysh formalism. More specifically, 
our main results are as follows: 

• We demonstrate that finite temperature strongly 
affects the screening properties of graphene, which, 
in turn, influences the inelastic scattering rates as 
compared to the zero-temperature case. 

• Focussing on the high-temperature regime, we cal- 
culate the quantum scattering rate which is rele- 
vant for dephasing of interference processes. We 
identify an hierarchy of regimes arising due to the 
interplay of a plasmon enhancement of the scatter- 
ing and finite-temperature screening of the interac- 
tion. 

• We further discuss the energy relaxation and trans- 
port scattering rate in graphene. Our results com- 
plement the kinetic-equation and hydrodynamic 
approache d 56 ' 57 ' 59 for the collision-limited conduc- 
tivity of clean graphene. 

The paper is organized as follows. In section |H] we 
define the model and develop the Keldysh diagrammatic 
formalism for treating the problem of Coulomb interac- 
tion in graphene. Section [Hi] is devoted to the analysis of 
polarization operator of graphene. We first review and 
discuss the results for the zero temperature case. Then 
we turn to the case of finite temperature and discuss 
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the properties of the polarization operator. In particu- 
lar, we compare the approximate form of the dynamically 
screened interaction propagator with the exact numeri- 
cal results. In Sec. IIVI we calculate the inelastic scatter- 
ing rates in the random phase approximation. Here we 
use the finite-temperature polarization operator obtained 
in Sec. IIIII to treat the problem analytically. We dis- 
cuss the asymptotics of the inelastic rates and show that 
Dirac fermions show no Fermi liquid behavior. Further, 
in Sec. [V] we discuss the collision-limited conductivity 
obtained within the diagrammatic approach. In Sec. IVII 
we conclude and summarize the main results of this pa- 
per. Technical details of the calculations are presented 
in three Appendices. 



II. THE MODEL 



A. Clean graphene with Coulomb interaction 



In this paper we consider clean graphene near the de- 
generacy point. The problem is described by the follow- 
ing Hamiltonian, which is a sum of the Dirac Hamiltonian 
Hq (describing the physics of non-interacting electrons in 
graphene at not too high energies) and the Coulomb in- 
teraction term, 



H = H + V=Y] [ d 2 ri>i(v)(~iv F a^)i> v {v) + \ V / d 2 n d 2 r 2 $t(n) &t,(r 2 ) —1. r fc„(r 2 )$„(ri 

„ J ^ — ,J £1*1 — r 2 



(1) 



Here e is the dielectric constant. The spinors \& have 
two components in the sublattice space, at are Pauli ma- 
trices operating in this space. The indices v, v' label N 
independent degrees of freedom (in graphene N = 4 ac- 
counts for spin and valleys degeneracy): the Coulomb 
interaction is invariant with respect to any rotations in 
the corresponding space. We set ft = 1. We focus on the 
case of undoped graphene and set the chemical potential 
(counted from the Dirac point) to zero, /i = 0. 

The retarded (advanced) Green's function of the non- 
interacting Hamiltonian Ho (the bare Green's function) 
in the energy-momentum space has the form 



Go ( e >P) 



el + v F cr ■ p 

(e ± ioy 



v 2 F p 2 



(2) 



It is convenient to introduce the projection operators 
that distinguish between the two chiral states: 



-P±(P) = 



t ± a ■ n 



(3) 



where n p = p/p is the unit vector in the direction of mo- 
mentum. With the help of Eq. ^ the matrix Green's 
function, Eq. @, can be decomposed into the superpo- 
sition of the two Green's functions corresponding to the 
states with + and — chiralities: 



Gf' A (e,p) = V + (p)G 1 f(e,p)+V-(p)Gt A (e, P ) , (4) 



and G A ±(e, p) = [Gffe.(e, p)] *. For the later purposes we 
will need the quasiparticle spectral weight 

A (e,p) = ±[GK(e,p)-G$(e,p)] = 

- ^ (el + v F (T ■ p)[S(e - v F p) + S(e + v F p)} . (6) 

Using the projection operators, we decompose the spec- 
tral weight as follows: 

■4o(e,p) = V + {p)A 0+ (e,p)+V-{p)A -(e,p), (7) 
A±(e,p) = -TrS{eTv F p). (8) 

It is worth noticing that for Dirac particles the spectral 
weight Ao(e, p) is not given by the imaginary part of the 
Green's function Gg : (e,p), because the latter contains 
the Pauli matrix a y . However, within each chirality the 
conventional relation holds: Aq± (e, p) = Im Gq± (e, p) . 

Next, we introduce the coupling constant for Coulomb 
interaction in graphene 



a a 



e 

ev F ' 



(9) 



which is similar to the fine structure constant 

e 2 1 



a = — 



137 



but is c/ev F times larger. The bare propagator of the 
Coulomb interaction in (fTJ reads (in momentum space): 



A)(q) 



2ira q v F 



(10) 



where 



Go ± («>P) 



1 



e + iO =F v F p 



(5) 



Throughout the paper we assume a g -C 1. This assump- 
tion is favored by recent experiments 27 ' 34 which sug- 
gested that the effective interaction constant in graphene 
is rather small at experimentally relevant temperatures. 
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B. Keldysh formalism 

Although in this work we discuss only equilibrium 
physics, we use the Keldysh formalis m 63 ' 64 in order to 
have a basis which can be generalized to the nonequi- 
librium situation. In particular, this Keldysh formalism 
will be used elsewhere for deriving kinetic equations for 
clean and disordered graphene. 

The bare Green's function is now a matrix in the 
Keldysh space 



(11) 



Here the Keldysh component Gq at the equilibrium 
reads 

C#(e,p)=/(e)[G?(e,p)-<#( e ,p)] 

= 2i/(e)^ (e,p), (12) 

where the fermionic thermal factor is given by /(e) = 
tanh(e/2T). 

The full Keldysh Green's function for Eq. ([T]) is ex- 
pressed through Go 



G = (G -E) 1 



(13) 



where S is the full self energy. In the lowest order in the 
fully dressed propagator of Coulomb interaction D the 
retarded self-energy is given by 



(14) 



where the symbol o denotes integration over all internal 
energies and momenta. 

In the equilibrium situation the Keldysh component 
D K of the interaction propagator satisfies 



D K {uj,q) = 2ig(Lu)lmD H (io,q) , 



(15) 



where g(to) = coth(ti;/2T) is the bosonic thermal factor. 
Using Eqs. (flip. (fT%|) . and (|T2]), one gets 



(16) 



-2i [(/ + g) Im D R ] o (G R - G J 



The retarded self energy is a matrix which contains the 
two terms: 



Y, R = Ef 1 + E* a ■ n v 



(17) 



The real parts of Ef and Ef give rise to the corrections 
to the energy and Fermi velocity, respectively: 



Se = -ReEf(e,p), 
Sv F = Re£?(e,p)/p. 



(18) 
(19) 



In order to fix the energy unrenormalized, we introduced 
the Z factor 



Z(e,p) = l 



ReEf(e,p) 



(20) 



so that e — ReEf (e, p) = Z(e,p)e. The renormalized 
Fermi velocity takes then the form 

ReE«(6,p)" 



Vp = v F Z(e,p) 



1 



(21) 



Using Eqs. ([2^)1 and ([2"Tjl , we write the retarded Green's 
function as 



G fl (e,p) 



= Z 



(e-iZImEf)l 



iZImE« 
P 



a p 



where 



ImEf = -~lmD R (f + g)oTrA, 



ImE. 



--ImD R (f + g)oTrA(T-n TJ . 



(22) 

(23) 
(24) 



Note that in this representation the velocity acquires a 
non-zero imaginary part. 

The most singular terms in the renormalized velocity 
v F , coupling constant a* = e 2 /ev* F , and the Z- factor can 
be summed up by means of the renormalization group 
approach 3 ^. For the case of weak interaction a g <§; 1 the 
solution of one-loop RG equations has the form: 



v F (e) 
Z(e) 



= Vp 1 



a g A 
-r- In — 



■ In ■ 



exp 



a 2 In- 



3tt 1 + 2l In : 



(25) 
(26) 

(27) 



Here A is the ultraviolet energy cutoff (bandwidth) and 
the on-shell relation between e and p is assumed. At finite 
temperature T the renormalization stops at max[e, T]. 
Therefore, for energies below T (which will be in the focus 
below) the renormalized velocity and the Z-factor are 
independent of energy. Since for a g <C 1 the corrections 
to unity in the renormalized Z-factor are parametrically 
small, in what follows we set Z = 1. 

Similarly to the bare Green's function, Eq. (Q| , the full 
Green's function, Eq. (|2"2"|) , can be represented as a sum 
of the two terms corresponding to ±-chiralities: 

G R (e,p) =V+(p)G R (e,p) +V-(p)G R (e 7 p) 7 (28) 



where 



with 



G£(e )P ) 



e T v F p - 



iZImE^ 



(29) 



ImEf = ImEf ± ImEf 



lmD R (f + g)oTr[AV±}. (30) 
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Clearly, the bare Coulomb interaction (whose propa- 
gator is purely real) does not yield the imaginary part of 
the self-energy, so that one has to take into account the 
retardation effects. In the random phase approximation 
(RPA) , the screened Coulomb interaction takes the form 



£>RPA(w,q) = 



D (q) 



1 + D (q) MI*(w,q) 



(31) 



where N is the number of flavors. The dynamical screen- 
ing in Eq. pip is expressed through the bare polarization 
operator 

n R = l -Tr(GXo G « + G«oG$). (32) 

In Sec. IIIII below we study II at zero and finite tempera- 
ture. 



III. POLARIZATION OPERATOR 

In the present section we discuss the properties of the 
polarization operator. This is of primary importance for 
understanding the screening of the electron-electron in- 
teraction and thus the physics of interaction-induced phe- 
nomena. We find it instructive to start with analyzing 
the zero temperature result and describing the processes 
relevant for the polarization operator. Then we turn to 
the case of finite temperature, which is our main interest 
in the paper. Finally, we analyze the consequences for 
the RPA-screened interaction. In the end of this section 
we comment on the applicability of the RPA in graphene. 



A. Polarization Operator at Zero Temperature and 
RPA Interaction Propagator 

The polarization operator in the energy-momentum 
representation reads 

n^,q)^-/|^/|/(e)Tr Mo ( e ,p) 

x [G£(e + W ,p + q) + G^(e-w,p-q)]} , (33) 

where G : are bare Green functions (J2J) . The momen- 
tum integrals that appear in this expression can be conve- 
niently evaluated using elliptic coordinates as described 
in Appendix [A"l 

In the zero temperature limit, Eq. (IA1[) simplifies and 
leads to the well-known expression for the imaginary part 
of the polarization operator^ 



-,/;, ._ !. Re . 9 2 signw 



Im II K (w,q) 



16 



2 _2 



(34) 



v F q 



One can see that the imaginary part of the polarization 
operator is non- vanishing only if uj > vpq and shows a 
divergence at the "light cone" uj — VFq- 



To understand the vanishing of the imaginary part of 
the polarization operator at u < vpq, it is instructive 
to analyz e) 54 ' 59 the kinematic restrictions for elementary 
processes. It is easy to see that an on-shell electron-hole 
pair can be created if (left panel of Fig. [T]) 



uj > VFq ■ 



(35) 



This should be contrasted to the energy and momentum 
conservation of on-shell electron-electron scattering pro- 
cesses (right panel of Fig. [1) 

uj 2 + 2v F p lP2 (l - cos(<[pi, p 2 ])) = v 2 F q 2 , (36) 

which implies the condition 

uj < vpq. (37) 

Thus, if we restrict our consideration to on-shell par- 
ticles, electron-electron scattering processes only create 
electron-hole pairs under the condition uj = VFq, when 
Imll diverges. As follows from Eq. (|36|) . scattering pro- 
cesses that satisfy uj — VFq correspond to forward scat- 
tering with <[pi,P2] = 0. After the RPA resummation, 
the interaction propagator takes the form 



2ttvfo.„ 



TTVFOiqN 



(38) 
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V^ 2 



v'pq 2 



It is seen that the divergence at uj = VFq is now elimi- 
nated; moreover, the imaginary part of the propagator is 



e-h 



e-e 





FIG. 1: Schematics of electron- hole (e-h) creation (left panel) 
vs electron-electron (e-e) scattering (right panel) near the 
Dirac point for /i = 0. The hole-hole scattering processes are 
analogous to those shown in the right panel. We set vp = 1 
for brevity. In both panels the possible transferred momenta 
satisfy |o — 6| < |q| < a + b (only the extreme cases of minimal 
and maximal q are shown). For e-h creation u> — a + b, so that 
w > qvF, whereas for e-e scattering the kinematic restrictions 
yield uj — a — b < qvF- 
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FIG. 2: Interaction-induced dispersion correction near the 
Dirac cone. The electron- hole creation [oj > vp(q)q] and 
electron-electron scattering [u < vf(p) q] regions do not over- 
lap, so that the zero-temperature inelastic scattering rate is 
zero. 



zero at the light cone. Therefore, the forward scattering 
divergence that arises on the Golden Rule (GR) level dis- 
appears within RPA, yielding a zero scattering rate on 
the RPA level. This is a manifestation of a highly sin- 
gular character of the zero-temperature problem where 
RPA may be insufficient. We will see below that at fi- 
nite temperature the thermal broadening regularizes the 
problem, so that scattering rates can be evaluated within 
the RPA. 

In view of the singular character of the problem, at 
zero temperature the imaginary part of the self-energy is 
highly sensitive to changes of the electron dispersion 59 . 
The Fermi velocity depends logarithmically on the mo- 
mentum or energy^ 5 -, Eq. (1251) , due to the renormalization 
by Coulomb interaction. In Eq. (f34f we have neglected 
the momentum dependence of vf which leads to a sepa- 
ration of the two regions defined by Eqs. (1551) and (1571) . 
see Fig. [5J The electron- hole creation [ui > vp(q)q) and 
electron-electron scattering [ui < vp(p) q] regions are then 
separated due to the renormalization-induced nonlinear- 
ity of the electron dispersion, which leads to the vanishing 
of the zero-T scattering rate already at the GR level. 

However, at finite temperature the situation is essen- 
tially different. First, the conditions ([35]) and (|37| will 
be smeared by temperature. Second, for energy scales 
smaller than T the renormalization of Fermi velocity 
is cut off by temperature and hence, the linearity of 
the dispersion relations is restored: renormalization re- 
duces merely to a T-dependence of the Fermi velocity. 
Therefore, when discussing the finite-T physics of in- 
elastic scattering on scales < T we will disregard the 
renormalization-induced nonlinearity of the dispersion. 
Furthermore, below, whenever we will use the notation 
vf for e < T we will mean the renormalized value of the 
velocity v* F (T) and omit the asterisk for brevity. 



B. Polarization Operator at non-Zero Temperature 

We are now ready to calculate the polarization opera- 
tor at finite temperature. It turns out that the effect of 
finite temperature on the screening in graphene is much 
more pronounced than in conventional metals with fi- 
nite Fermi-surface and quadratic electronic dispersion. 
Indeed, the linearity of the spectrum of Dirac fermions 
gives rise to a strong (linear) energy dependence of the 
density of states in two dimensions, whereas for parabolic 
spectrum the density of states is constant. In the latter 
case, the polarization operator is essentially independent 
of temperature. This is not the case for Dirac particles. 

Physically, a finite temperature leads to population 
of electronic states in an energy range ~ T around the 
Fermi level. Let us consider the undoped graphene where 
the chemical potential lies at the Dirac point. The typi- 
cal density of states that participate in the screening of 
Coulomb interaction at finite temperature is now pro- 
portional to T. In more technical terms, at finite T 
the integration over the fermionic energy in the polariza- 
tion bubble essentially involves not only the distribution 
function [thermal factors /(e)] but also the Green's func- 
tions. This strongly changes the polarization operator at 
qvF,u < T. 

To simplify the notations at non-zero temperature, we 
introduce the dimensionless variables according to 



VFq 
2T 



2T 



(39) 



We use the general expression for the polarization oper- 
ator, Eqs. (|A1| and (|A2j) from Appendix lAl Considering 
four different regions shown in Fig. [3J we simplify these 
equations in each of the cases, which allows us to treat 




FIG. 3: (Color online) Imaginary part of the polarization op- 
erator in the frequency-momentum plane. Four regions of 
dimensionless variables Q and Q (see text) are indicated. 
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the problem analytically. • Region 3: Q <C 1 and |fi| < Q 

The four regions are defined by the following condi- 
tions- • Region 4: 1 < Q and |0| < Q 

• Region 1- Q <C 1 and Q < Ifil ^ s snown m Appendix El the leading-order terms with 

non- vanishing imaginary part form the following simpli- 

• Region 2: 1 C Q and Q < |f2| fied polarization operator: 



n fl (tt,Q) 



T 
775" 



1 - 



In 2 

7T 

In 2 " / In 2 

7T 



mi 



v^ 2 - Q 2 



i Q 2 taah(fi/2) 

8 V^-Q 2 



i Q signil 



8 ^-Q 2 
1 Q 2 



Qe-^sinhO 



8 ^Q2 - ft2 ^ ^/Q2 _ ft2 

I 



(re* 


;ion 


1) 


|n| 


> Q 


(re{ 


;ion 


3) 


|n| 


< Q 


(re{ 


,ion 


2) 


|n| 


> Q 


(re§ 


,ion 


4) 


|n| 


< Q 



Q< l 



> 1 



(40) 



In the case of large momenta, Q 3> 1, we recover the 
zero-temperature result. For Q -C 1, the polarization 
operator substantally differs from the zero-T expression. 
We will discuss this case in more detail in section IIII CI 
The separation between < Q and fl > Q is dictated by 
the non-analytical structure of the polarization operator, 
see Eq. (|4T)]) . 



C. RPA Interaction at non-Zero Temperature 



The real part of the polarization operator, Eq. (|40| . for 
Q -C 1 is determined by temperature and leads therefore 
to the screening of Coulomb interaction: 



HmD&, A (n,Q) = - 



D (Q) 



D (Q)NIL R (Q,Q) 

VpTT 



a g N ln2 



NT\n2 \Q 
which yields the screening length 

^scr — 



a g N In 2 



Vp 



2a g NT\n2 



(41) 



(42) 



Thus at finite temperature the system does screen the 
long range Coulomb interaction. Note that for Q <C 
a g N, the RPA propagator becomes independent of the 
interaction constant a a , 



D* VA (n,Q) 



1 



Q«a s Ar - n«(o,Q)' 



(43) 



as in conventional systems with Coulomb interaction. We 
will see, however, that the dominant contributions to re- 
laxation rates are determined by higher transferred mo- 
menta Q > ct g N, where the peculiarities of the finite-T 
screening in graphene are crucially important. 



In Region 1, the real part of the polarization operator 
in Eq. (|40[) may become negative, leading to emergence 
of plasmon excitations. The plasmon dispersion Q p (Q) 
is determined by the zero of 



1 + D NU R = 



a g NQ In 2 



(fi 2 - Q 2 ) 



2^3/2 



n- 



(a g N\n2 + Q)^Q 



a g N2 In 2 



yielding 



{a g N\n2 + Q)^Q 



2a g N\n2 



(44) 



(45) 



For Q — > 0, this simplifies to Q, P (Q) cx ^Ja g NQ. A 
non-zero imaginary part of the polarization operator (|40j) 
in the corresponding region implies that these plasmons 
have a finite lifetime. The decay rate of plasmon excita- 
tions is given by 



r„(fi,Q) 



(n 2 - Qr 



.3/2 



aJVQ In 2 



■lmD {Q)ImU R (n,Q), (46) 



which yields 

r p (Q, Q)| n= n p (Q) 



16 In 2' 



n 2 



< 



-Q 2 



16 In 2 



n=n„(Q) 



n=n p (Q) 



(47) 



Remarkably, Eq. (|4"T|) indicates a good quasiparticle be- 
havior for plasmons which is, as we will see in Sec. IIV A[ 
not true for electronic excitations. The situation is some- 
what similar to Luttinger liquid where plasmons are al- 
most perfect quasiparticles, whereas, from the point of 
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FIG. 4: (Color online) a) Imaginary part of RPA interaction. Left: using an analytic approximation, right: by numerical 
evaluation of Eqs. (lAljl and (|A2[) . Dashed lines indicate the plasmon dispersion (see Eq. ([44} ) . b) crosssections through the 
plots of a) from (fi, Q) = (1,0) to (fi, Q) = (0,1) (1) from (fi, Q) = (0.2,0) to (fi,Q) = (0,0.2) (2). Solid curves: analytic 
approximation; dashed curves: numerical integration. For all plots the interaction strength is a g N — 1/3. 



view of fermionic excitations, the Luttinger liquid repre- 
sents a paradigmatic example of a non-Fermi-liquid. 

Figure |4] demonstrates that Eq. (|40|) for the polariza- 
tion operator yields a remarkably good approximation for 
evaluation of the imaginary part of the RPA-interaction, 



Im Drpa 



NDl ImLT^ 



(1 + ND HjeU R ) +(ND lmU R ) 



(48) 



[see also Eq. (|B2|) in Appendix IB"] . In Fig. 2]d we see the 
plasmon peak in the RPA interaction propagator, which 
is strongly asymmetric and is suppressed around the light 
cone, so that ImZ?FiPA = exactly at 51 = Q. 

Let us now discuss the status of the RPA in graphene. 
Above we have introduced the large number of indepen- 
dent flavors N 1. At zero temperature, this is the 
only parameter which justifies the RPA summation in the 
problem of interacting Dirac fermions at /i = 0. Indeed, 
in view of the absence of the screening, the non-RPA di- 
agrams are parametrically the same as those included in 
the RPA series for N ~ 1. Furthermore, as discussed in 
the end of Sec. 1111 Al the renormalization-induced cur- 
vature of the spectrum (which is also beyond the RPA) 
may dramatically affect the results obtained within the 
RPA. 

However at finite T, the T-induced screening of the 
interaction, Eq. (|4"Tj) . restores the validity of the RPA 
for q < T/v F (Q < 1) even for N ~ 1. Indeed, the 
1/q-singularity of the long-range Coulomb interaction is 



not compensated in the denominator of Eq. (|48|) because 
the polarization operator at finite T is no longer linear- 
in-g. The situation becomes similar to that in conven- 
tional metals with a finite Fermi-surface, where the RPA 
is justified for q <C kp. In graphene the role of kp is 
played by T/vp, which in effect establishes an analog 
of a finite Fermi-surface. Therefore, for q <C T/vp the 
RPA does sum up the most singular interaction-induced 
terms: all other terms are non-singular because of the 
screening. This means that all the observables that are 
dominated by the collisions with the momentum transfer 
smaller than T jvp can be calculated (even with the cor- 
rect numerical prefactors) within the RPA. The RPA re- 
sult for those observables that are dominated by qvp ~ T 
is parametrically correct, but the value of the pref actor 
can not be found using the RPA. 

Below we employ the finitc-T RPA for calculation of 
various scattering rates in graphene. For the sake of gen- 
erality, we keep N as a parameter. In what follows we 
focus on the case a g N <C 1, but whenever the rate under 
the consideration is dominated by qvp < T, the condi- 
tion N 3> 1 can be removed so that we are allowed to 
use the RPA for N ~ 1. 



IV. SCATTERING RATES 

In this section we calculate various inelastic scatter- 
ing rates in clean graphene at finite temperature. More 
specifically, we focus on the quantum scattering rate (and 
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dephasing) and the energy relaxation rate induced by the 
RPA-screened Coulomb interaction. In Sec. IV Al below, 
we will also calculate the transport scattering rate due to 
inelastic collisions. The quantum scattering rate deter- 
mines the lifetime of quantum states (plane waves) and is 
related to the dephasing rate. The energy relaxation rate 
governs the relaxation of the quasiparticle distribution 
function. Finally, the transport scattering rate describes 
the influence of the inelastic scattering on transport phe- 
nomena. 

Although the origin of all these rates is the same — 
the inelastic electron-electron collisions, these rates, may 
strongly differ from each other. For instance, this is ex- 
actly what happens in diffusive metal o 33 ' 65 because of 
the infrared-singular collision kernel. Another prominent 
example of a non-trivial behavior of relaxation rates re- 
lated to the infrared singularities is a Luttinger liquid 
(disordered or clean) £§r£! On the other hand, within 
the Fermi-liquid theory of clean metals all the inelas- 
tic scattering rates behave in the same way, since the 
characteristic frequency /momenta transfer in the course 
of electron-electron collisions is determined by tempera- 
ture. The goal of this section is to understand whether 
the situation in clean graphene is similar to the Fermi- 
liquid or not. 

We have already mentioned in the introduction (Sec. 
HJ the previous works on the inelastic quantum scattering 
rate in graphen o 35 i 55 . These works addressed the scatter- 
ing rate at the GR level at zero temperature and obtained 
the Fermi-liquid- type result r" 1 ~ a?iVe <C e. Naively, 



C^N 2 T 



q,NT 



FIG. 5: Characteristic energy scales separating domains of 
distinct behavior of the rates. 



A. Quantum scattering rate in graphene 

One of the main manifestations of the inelastic scat- 
tering in electronic systems is the interaction-induced de- 
phasing. In order to analyze the effects of dephasing in 
graphene, we will follow the route suggested by earlier 
works on 2D (in particular, diffusive) systems.— A natu- 
ral first step is to calculate the quantum scattering rate 
t~ 1 which is given by the imaginary part of the quasipar- 
ticle self-energy taken at the mass-shell. Indeed, for con- 
ventional metals with parabolic dispersion, in the high- 
temperature (ballistic) regime, the dephasing rate t^ 1 
to leading order is given by r" 1 , see Ref. [7(J In the dif- 
fusive regime, the scattering kernel acquires an infrared 
singularity leading to a divergent r" 1 at finite T. This 
is a manifestation of the fact that the single-particle self- 
energy is not a gauge- invariant object; no divergencies 
occur in observable (gauge-invariant) quantities, such as, 
e.g., dephasing rate. However, even when the quantum 
scattering rate diverges, the calculation of it turns out 
to be instructive: a parametrically correct result for the 
dephasing rate can be obtained from the expression for 
t~ 1 supplemented with an appropriate infrared cutoff. It 



one could think that at finite temperature this consider- is tllus useful to be S in with analyzing r q 



ation would lead to r q 1 



-NT. However, since the 



GR result is completely determined by the "mass shell" 
(cj = vfq) (see discussion in Sec lIII A|) one concludes that 
the RPA-resummation (which kills the on-shell interac- 
tion) would yield t~ (e) = for T — and arbitrary e. 
Therefore, the finite- T expectation based on the GR is 
also doubtful. 

How the quasiparticle broadening behaves in the ex- 
perimentally relevant case of finite temperature is thus 
by far not obvious. As discussed in Sec lIII Bl the finite 
density of states at non-zero T leads to screening of the 
Coulomb interaction, thus justifying the use of the RPA 
which then sums up the most singular contributions of 
the interaction for small momenta, qvF < T. In this 
section we find the behavior of inelastic rates at finite 
temperature within the RPA. As we have already seen in 
section llllBI the behavior of the polarization operator at 
finite temperature is highly non-trivial. This leads to a 
rather rich behavior of the scattering rates. 

It turns out that, in addition to the temperature scale, 
two more characteristic scales appear which are relevant 
for relaxation rates: a 2 g N 2 T and a g NT. In this section 
we will distinguish between the four regimes (TIV) as 
shown in Fig. [SJ The contributions to each of the scat- 
tering rates from different regions (1-4, see Fig. [3]) of the 
momentum-frequency plane are calculated separately. 



1. Quantum scattering rate: definitions 

The peculiarity of graphene is that the self-energy is 
a matrix in the sublattice space, which has a non-zero 
imaginary part, Im E„ [Eq. (|24[)]. in the off-diagonal com- 
ponents. Therefore, the definition of the quantum scat- 
tering time in graphene is actually not unique. Indeed, 
one can associate with the quantum scattering rate the 
on-shell value of the imaginary correction to the energy 
in the full Green's function, ImS e [Eq. ([23]) ]. similarly to 
the conventional Fermi-liquid theory: 



1 



2r q (e) 



(e) ImS e (i; F p,p)| p= 



z/v F 



-6(-e)lmZ e (-v F p,p)\ p= _ e/vF (49) 



Here at e > we have taken the self-energy at the "+" 
mass-shell corresponding to the positive energies, e = 
vfP, and at e < on the "— " mass-shell. Clearly, 



1 



1 



i(e) 7q(-e) 



(50) 



for undoped graphene (// = 0) because of the particle- 
hole symmetry. Within the RPA the explicit expression 
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for the imaginary part of the total self-energy S e taken 
at e = v F p reads: 

7^ [g(u) + f(v FP -uj)} 



d 2 q 



ImD^u.q) 



x [S(v F p — bj — v F \p - q|) + S(v F p - uo + v F \p - q|)] . 

(51) 

Alternatively, one can introduce the lifetime of the 
+ and — chiral states through the corresponding self- 
energies, Eq. ([30|) : 



1 



2r±(e) 



= - ImE±(e,p)| p=|e|M 



which yields 



-(e) 



-(-*) 



(52) 



(53) 



Using Eqs. (|30]) and (|6|), we get for the self-energy of 
electrons (+ chirality) 

7T f d,U! 



d 2 g 



1 - 



(2tt)2 

i_ p(p - q) 
p|p - q| 
p(p - q) 
p|p-q| 



lmD H (uj,q) 
5(v F p — uj — v F \p — q 
8(v F p - ui + v F \p - q|) 



(54) 



for the "own" mass-shell, and 



ImS + (-w f p,p) = -| J ^-[g(uj) + f(-v F p-uj)} 



d 2 q 



1 - 



(2tt)2 

p(p - q) 
p|p-q| 
p(p - q) 
p|p-q| 



lmD R (u,q) 
5(—v F p — u) — v F \p — q|) 
S(-v F p -uj + v F \p - q|) 



(55) 



for the "wrong" (hole) mass-shell. 

The main formal difference between the two relaxation 
rates, Eq. (l49l) and Eq. (1521) . which are related by 



1 



1 



1 



r q (e) 2r + (e) ' 2r_(e)' (56) 

is in the appearance of Dirac factors 1 ± cos 9 in Eqs. (|54|) 
and (|5"5")l where = arccos(pp'/pp') is the scattering an- 
gle between the incoming momentum p and the momen- 
tum p' = p — q after scattering. These factors are related 



to the additional Berry phase in the problem of Dirac 
particles, which arises due to the overlap of Bloch func- 
tions and, in particular, forbids the backscattering within 
the same chirality and valley. Note that for well-defined 
quasiparticles (i.e., in a Fermi-liquid situation) the self- 
energy at a "wrong" mass shell would never be relevant. 
However, if the quasiparticle broadening is larger than 
the characteristic energy, this is no longer the case, so 
that Eq. (|55"|) may then contribute to the observables. 

In Eq. (f54"|) , the term with 1 + cos 9 corresponds to the 
electron-electron scattering (right panel of Fig. [T|) which 
is determined by the contribution of region 3 in Fig. [3] 
The term with 1 — cos 9 is due to electron- hole scattering 
(electron-hole pairs, left panel of Fig. [1]) accompanied 
by the excitation of plasmons and is determined by the 
contribution of region 1 in Fig. [3] The latter contribution 
is suppressed for the forward scattering 9 = because 
of the Dirac factor. Furthermore, at zero temperature 
only the electron-electron processes are allowed by the 
kinematic restrictions. 

At finite T, however, the situation is different for the 
low-energy domains I and II (e < a g NT), where the con- 
tributions of the electron-electron and electron-hole scat- 
tering to the inelastic quantum scattering rates are of the 
same order. This means that the low-energy electron- 
type and hole-type quasiparticles are strongly correlated 
by the mutual inelastic scattering, whereas at higher en- 
ergies (e > a g NT, corresponding to domains III and IV) 
the electronic and hole subsystems are only weakly cou- 
pled with each other, in agreement with Ref. [59l 

As we will see below, depending on the energy range, 
Fig. [5l both r" 1 and r^ 1 may be larger or smaller than 
energy. When the quasiparticle 's broadening is small 
( "Fermi-liquid regime" ) , the two rates coincide since the 
inelastic scattering is dominated by small scattering an- 
gles 9 and the Dirac factors reduce to and 1. At the 
lowest energies (domains I and II), both r^ 1 and r^ 1 ex- 
ceed the energy which makes the notion of mass shell not 
well-defined. In particular, due to the inelastic broaden- 
ing, the electronic excitation (+ chirality) has tails at 
negative energies which overlap with the hole mass shell. 

Furthermore, as shown in Sec. IIV C[ the character- 
istic rate of the energy relaxation (energy mixing due 
to the diffusion over energy) is in the "non-Fermi-liquid 
regime" of the same order as the quantum scattering rate, 
so that the electronic excitations constantly explore the 
hole mass shell and vice versa. In this situation, the 
broadening of quasiparticles of a given chirality is in ef- 
fect described by the total scattering rate r" 1 , according 
to Eq. (|5o]) . rather than Eq. (f52"|) . Therefore, in what fol- 
lows, we will mostly focus on the Fermi-liquid-type total 
rate r" 1 , formally defined in Eq. (|4*^| . 



2. Quantum scattering rates: results and discussion 

The evaluation of integrals involved in the calculation 
of the total quantum scattering rate r" 1 is outlined in 
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FIG. 6: (Color online) Quantum scattering rate for a g N — 4x 
1CP 3 (double logarithmic scale). Dots: exact values obtained 



by numerical evaluation; solid lines: 

Eq. ung. 



analytical asymptotics, 
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II 


III 


IV 


qVF ~ cj 


T 


(a g NT) 2 /e 


Q 9 ivr 


a 9 JVT 


Regions 


1,2,3,4 


1,3 


3 


3 



TABLE 1: Momentum/frequency scales and regions of the 
Q-Q plane that dominate the quantum scattering rate t^ 1 in 
different domains (I, II, III, and IV) of energy e. 



Appendix [B] The result depends on the energy range 
(see Fig. 0): 



T U 



Cl N\ T 



e<a 2 g N 2 T, I 



2r q (e) 



c 2 a g T, a 2 N 2 T < e < a g NT, II 



c 3 a g T, 



a g NT<^e. 



Ill, IV 

(57) 

Here c\ ~ 1, C2 = 37r/2, and c 3 = 2/3(2) /n are the numer- 
ical coefficients of order unity, is the Dirichlet beta 
function and (3(2) is the Catalan constant. Note that for 
energies in Regimes II, III, and IV, the number of inde- 
pendent flavors N drops out from the expression for r q . 
The comparison of the asymptotic expressions (f57j) with 
the exact numerical evaluations is shown in Fig. [6]). 

The obtained rates are dominated by different values 
of momenta transferred during the collision. Further- 
more, depending on the energy range the main contri- 
bution may come from different regions of the Q vs. f2 
plane (Fig. [3]) as shown in Table [TJ 

In Regime I the result is determined by momenta of or- 
der temperature and therefore all four regions contribute 
in the same way: r^ 1 ~ (eT) 1 / 2 /N . All scattering an- 
gles 9 contribute to the result in Regime I. In order to 
evaluate the numerical prefactor c\ one needs the knowl- 
edge of the screened interaction in the crossover around 
qvF ~ T, which is beyond the analytic approximations 



r q 1 is domi- 



for the polarization operator used above. For energies 
in Regime II the dominant contributions come from re- 
gions 3 (electron-electron scattering) and 1 (electron-hole 
scattering); in both of them all scattering angles 9 con- 
tribute. In Regimes III and IV the main contribution 
to the quantum scattering rate stems from the region 
3 (electron-electron scattering) and is dominated by the 
forward scattering (9 < a g NT/e <1). 

An important feature of the quantum scattering rate 
is its non-monotonic energy dependence, see Eq. (|57|) 
and Fig. |5J with increasing energy the quasiparticle 
broadening first grows in Regime I, has a maximum at 
e ~ a 2 N 2 T, then decreases in Regime II, and finally be- 
comes energy-independent. The maximum of r^ 1 occurs 
due to the resonant emission/absorption of plasmonic ex- 
citations. We will see below that the non-monotonicity 
of the energy dependence is related to the peculiar prop- 
erties of the dynamical screening in graphene and is a 
characteristic feature of all inelastic scattering rates in 
graphene. 

In Fig. [7] we show the results for the quantum scatter- 
ing rate of electrons, r7 , defined in Eq. (|52[) . One sees 
that Fig. [7] represents an "unfolding" of Fig. [5] into the 
contributions of "own" (e > 0) and "wrong" (e < 0) mass 
shells, according to Eq. (|5"o| . The hole mass shell with 
e < is probed by electrons due to the strong quasipar- 
ticle broadening. At positive energies in Regimes III and 
IV (e > a g NT), the total scattering rate r~ 

nated by the electronic scattering rate r^I 1 . 

In contrast to the naive expectation (t^ 1 ~ a 2 NT) 
based on the Fermi GR result, the quantum scattering 
rate calculated within the RPA is proportional to a g (and 
does not depend on N) for e 3> a 2 N 2 T and is indepen- 
dent of a g at the lowest energies. This enhancement of 
the inelastic scattering for a g N <C 1 is a result of pe- 
culiar screening properties of graphene at finite T which 
leads to a nonanalytic behavior of the rate as a function 
of the natural four-fermion coupling constant a 2 . This 
behavior bears a certain similarity to the behavior of the 
quantum scattering rate in a spinful Luttinger liquidi 66 i 

Importantly, the formal calculation of the total scat- 
tering rate t^ 1 in the lowest order in the RPA inter- 
action propagator yields T~ 1 (e) > |e| for regimes I and 
II (e < a g NT). This means that the quantum scatter- 
ing rate is ill-defined in these low-energy domains. This 
also signifies that for the calculation of observables in 
these energy domains the higher-order terms in the RPA 
screened interaction may become important. Neverthe- 
less, the above lowest-order calculation in domains I and 
II is useful as it provides the characteristic value for the 
intensity of inelastic processes. 

Since the RPA quasiparticle broadening in domain I 
may overlap with the higher-energy domains, where the 
quantum scattering rate is of the order of a g T, one 
can speculate that the characteristic strength of inelastic 
scattering is given by T q " 1 (e) ~ a g T also at low energies. 
Of course, the calculation based on the lowest-order RPA 
diagrams for the self-energy is then insufficient. More- 
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FIG. 7: (Color online) Quantum scattering rate r + x for the + chirality (electrons) for a g N — 4 x 10 3 (double logarithmic 
scale). Dots: exact values obtained by numerical evaluation; solid lines: analytical asymptotics (see Appendix IB]! . 



over, the typical observables at finite temperatures are 
dominated by e ~ T (the border between domains III and 
IV), where the quasiclassical broadening is smaller than 
energy and the above calculation is justified. The low- 
energy inelastic relaxation may become relevant in the 
context of spectroscopy under strongly non-equilibrium 
conditions, for example, in problems related to tunnel- 
ing into a non-equilibrium state (cf. Ref. |69| ). In this 
situation the inelastic effects can be treated within the 
quantum-kinetic approach including non-RPA contribu- 
tions, similarly to one-dimensional problems; 69 



B. Dephasing rate at high temperatures 



as the resulting rate [see Tableland Eq. {57}], both are 
~ a g T for N ~ 1 (for N 3> 1, the characteristic mo- 
menta are much higher than l/r q ). In Regime I of lowest 
energies, the characteristic momentum transfer is much 
higher than the rate. Therefore, in all these regimes the 
infrared cutoff is redundant and the dephasing rate is 
given by the quantum scattering rate, Eq. lfST)) . similarly 
to the case of conventional metals in the ballistic regime 7 -^ 
(although, in contrast to the conventional case, the char- 
acteristic transferred frequencies are much smaller than 
temperature). Since the characteristic energies involved 
in the transport coefficients are of order of T, we con- 
clude that in the ballistic regime the interference effects 
are governed by 



Let us now discuss the dephasing rate. The result 
for the quantum scattering rate obtained for the clean 
graphene allows us to evaluate the dephasing rate rele- 
vant to weak (anti)localization in the ballistic regime^ 
of high temperatures Tt^s 3> 1, where r^is is the elas- 
tic mean free time due to scattering off impurities. It 
is worth noting that this condition may coexist with the 
condition 3> Tdi s which allows long interfering paths 
in the weak-localization experiment. Indeed, the inelastic 
scattering is suppressed with decreasing a g (even though 
in a non-trivial way) so that one expects that can be 
made arbitrary long. 

The result for the quantum scattering rate, Eq. (1571) . 
remains intact (up to the change in prefactor) when we 
calculate it in a self consistent way as appropriate for 
estimating the dephasing rate: 



1 



a g T. 



(59) 



This prediction can be verified by transport experiments 
on graphene at sufficiently high temperatures (depending 
on the purity of the system, T > 10 — 10QK for a typical 
setup with graphene deposited on a insulating substrate 
and T > 1 — 10K for suspended graphene flakes). At 
lower temperatures Tr^is <C 1 corresponding to the "dif- 
fusive regime" with respect to interaction, one expect the 
conventiona l 33 ' 71 diffusive result for the dephasing rate 



1 T 
In g, 

t<p g 



(60) 



where g is the dimensionless Dirac-point conductance. 
At the Dirac point g is close to unity and hence t7 1 ~ T. 



dQ. 



(58) 



This happens due to the fact that the characteristic mo- 
menta dominating the integrals for the quantum scatter- 
ing rate in Regimes II, III, and IV is of the same order 



C. Energy relaxation rate 

Let us now discuss the energy relaxation time in clean 
graphene. In Sec. IIV Al we have seen that the typical mo- 
mentum or energy transfer during the electron-electron 
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collision is much smaller than temperature. In this sit- 
uation, the energy relaxation occurs through multiple 
scattering processes which can be viewed as diffusion in 
energy space (see, e.g., Ref. 111). The characteristic en- 
ergy relaxation rate is given by the diffusion coefficient 
of this problem, which amounts to introducing the factor 
JCb = uj 2 /T 2 , see Appendix IB 11 into the collision ker- 
nel. More rigorous calculation of the energy relaxation 
or equilibritation rates can be done using the language of 
kinetic equation; here we only estimate the typical rate 
t^ 1 within the energy diffusion picture. 

Once the quantum scattering rate is obtained, the cal- 
culation of the energy relaxation rate can be done using 
the same steps as described in the previous section. Tech- 
nically, the integrals in Eq. ([5Tjl are only slightly changed, 
which leads, however, to a substantial difference between 
the two rates. The detailed calculation can be found in 
appendix [B] here we present only the result: 



T FT 
N V T' 



T E (e) 



e/T_ 
a 2 g N 2 



9 /e\3/2 
< NT (t) In 



a g N 



I 



II, III 



IV 



(61) 



Again the obtained rates are dominated by different mo- 
mentum scales and by contributions of different regions 
as shown in table [2] Due to the factor /Ce all contri- 
butions except in regime IV are now determined by mo- 
menta of order temperature, which does not allow us to 
find the numerical value of the prefactors analytically. 

Furthermore, in fact, the above calculation based on 
the energy diffusion is not justified for "hot electrons" 
with high energies, e > T . Indeed, within our consider- 
ation, the characteristic energy transfer dominating the 
energy relaxation in Regime IV turns out to be much 
higher than T, in contrast to the original assumption. 
Therefore, the estimate Eq. (ItTTI) for Regime IV can not be 
trusted and another approach is needed for this Regime. 
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FIG. 8: (Color online) Energy relaxation rate for a g N — 
4x 10 -3 (double logarithmic scale), obtained from the energy- 
diffusion consideration. Dots: exact values obtained by 
numerical evaluation; solid lines: analytical asymptotics, 
Eq. (|61[) . Since the energy-diffusion approximation employed 
in the calculation is not justified in Regime IV, we do not 
present the results in this Regime for e»T. 





Vpff 



FIG. 9: Diagrams describing the first-order interaction cor- 
rection to the conductivity 



resonant condition correspond to e ~ a 2 g N 2 T), whereas 
at high energies the energy relaxation is stronger because 
of the large phase space available for inelastic processes. 
We remind the reader, however, that at e > T the above 
calculation based on the energy-diffusion approximation 
is not justified. In order to find the correct relaxation 
of the distribution functions in this Regime of hot elec- 
trons, one should solve the corresponding kinetic equa- 
tion which can be reduced to the Fokker-Planck equation. 
This will be done elsewhere. 



TRANSPORT RATE AND CONDUCTIVITY 



TABLE 2: Momentum/frequency scales and regions of the 
Q-Q plane that dominate the energy relaxation rate r^ 1 in 
different domains (I, II, III, and IV) of energy e. 



The numerical results for the energy relaxation rate 
are shown in Fig. [5] together with the analytical asymp- 
totics, Eq. (joTj) . One sees that the energy relaxation has 
a minimum at e ~ T, where we recover, up to the log- 
arithmic factor 



lna s |, the GR result t e 



a 2 g NT. At 

lower energies the inelastic scattering is enhanced due 
to the resonance in the RPA interaction propagator (the 



Transport scattering rate due to inelastic 
collisions 



In this section we calculate the transport relaxation 
rate due to inelastic collisions. The expression for the cor- 
responding kernel of the self-energy can be deduced from 
the interaction-induced correction to the conductivity, as 
described below. We are interested in the linear-response 
dc conductivity. The leading order perturbative correc- 
tion to the conductivity due to Coulomb interaction is 
given by the two diagrams shown in Fig. [9l 

In the absence of interaction, the conductivity of a 
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clean graphene diverges at finite temperature. There- 
fore, we have regularized the diagrams by introducing 
a small broadening (<5) which mimics the finite lifetime 
due to weak disorder. The general analytic expression of 



the lowest-order interaction correction to the conductiv- 
ity within the Keldysh formalism can be found in Ref. 
l7ll . This correction can be split into two parts (Fig. [9^ 
and Fig. respectively): da = 8a s + 6a v , where 



1 f d 2 p! f dei f d 2 p 2 f de 2 dfi 



Tr \j [G R (l)-G A (l)]j f} 



6<J W At J (2ir) 2 J 2tt J (2tt) 2 J 2ir de 

x (/ 2 + g) [D R - D A ] [G R (l) G R {2) G R (1) - G A (1) G A (2) G A (1)] 
+ f 2 \G R (l)[G R (2)D R -G A (2)D A ]G R (l)-G A (l)[G R (2)D R -G A (2)D A ]G A (l)^ j (62) 
is the self-energy contribution and 

„ 1 f d 2 Pl fde 1 f d 2 p 2 fde 2 df 1 r. f R A] 

x [G R {1) G R {2) - G A (1) G A {2)\ j p [G R {2) G R {1) - G A {2) G A {1)\ 
+ 2 Jp h[G R {l) - G A (1)} [D A G R {2) ]p G R (2) - D R G A {2) j G A (2)} jp[G R (l) - G A (l)} } (63) 

is due to the vertex correction. Here jp = evF&p is the current operator in graphene and we used the short-hand 
notations for the arguments of Green's functions: 1 = pi,ei and 2 = p2,£2- We first trace out the sublattice 
structure and then regularize the divergent integrals by 8. The most divergent interaction-induced correction to the 
dc conductivity is then proportional to 5~ 2 : 



2, ,2 



5a 



e v 



AS 2 



dei df f de 2 f d 2 pi f d 2 p 2 



2tt dei J 2tt J (2tt) 2 J (2tt) 



7r[<5(ei - v F pi) + S(ei + v F px)] signei 



1 - 



P1P2 

P1P2 



Pi ' P2 

PlP2 



K s (ei,ei - e 2 ,Pi,Pi - p 2 )- (64) 



Transport factor 



Dirac facto 



Here is the integral kernel of the self-energy: 

ImE«(e,p) = J-^L |^K s ( e , w ,p,q) (65) 

The correction, Eq. is determined by the inelastic 
electron-electron scattering. In the diagrams of the lead- 
ing order in bare interaction shown in Fig. [SJ we have 
Kg cx ImDo = 0, so that the inelastic corrections are 
zero. What remains in the first-order conductivity cor- 
rection, Fig. |9l are the contributions responsible for the 
renormalization of the Fermi velocity in graphene com- 
ing from the real part of the self-energy; note that these 
contributions are less singular in <5 _1 . 

Thus, we have to consider higher order corrections. In 
Fig. [10] one can see all classes of second order skeleton 
diagrams. There are also the second-order diagrams of 
the ladder type, see Fig. [IT] Such diagrams contribute 



only to the renormalization of Dp, similarly to the dia- 
grams in Fig. [9] and their contribution has been already 
included into the calculation simply by the replacement 
v F -> v F (T). 

In the large N approximation diagrams in Fig. llOt i. 
[TDfe . andlTUf dominate. DiagramlTOf is known as Coulomb 
drag diagram and yields a zero contribution at the Dirac 
point due to electron-hole symmetry (see Refs. [73lf74l and 
Appendix [Cjl . Thus we are left with diagrams [TDH and 
[TOe which correspond to the two diagrams shown in Fig. 
[9] with the second-order correction to the interaction in- 
stead of the bare one. For these two diagrams Eq. (|64|) 
still holds with Ks oc ImDi. This suggests to replace 
the bare interaction lines shown in Fig. [9] by the RPA in- 
teraction lines, which would correspond to Eq. (|64[) with 
Kj] cx IiiiDrpa. Note that even with the RPA-dressed 
interaction, Eq. (1641) still yields a divergent contribution 
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FIG. 10: Second-order skeleton diagrams for the conductivity 
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FIG. 11: Examples of second-order diagrams for the conduc- 
tivity that contribute to the renormalization of velocity. 



which we have regularized with S. 

Considering the two independent scattering processes 
with the transport scattering rates S and t^ 1 such that 
we find the Drude conductivity in the form 



' J dep(e) (- 



dn,F \ VpT t r(e) 
de J 



2[1 + T tr (e)<5]' 



(66) 



where njr(e) = [1 — /(e)]/2 is the thermal Fermi distri- 
bution function and p{e) is the thermodynamic density 
of states. Expanding this formula in (5r tr ) _1 , we get the 
"interaction- induced" correction: 



5a 



,2„,2 



AS 2 



de 



df p(e) 



de T tr (e) 



(67) 



Comparing Eq. (I67p with Eq. (|64[) we can iden- 
tify the interaction-induced transport scattering rate in 
graphene. The transport scattering rate obtained from 
the conductivity correction corresponds to the kernel /Ct r 
defined in Appendix IB 11 

Furthermore, the connection between the kernels in 
the transport scattering rate and the quantum scatter- 
ing rate also follows from Eqs. (|6"51 . and ([BT)) . One 
can see from Eq. that in the transport scattering 
rate not only the contribution of forward scattering pro- 
cesses is suppressed as in conventional systems but also 
the contribution of the backward scattering. The latter 
suppression is due to the Berry phase of n in graphene. 

In the above derivation we have assumed S 3> T t ~ 
which allowed us to extract T t ~ from the expansion of 
the conductivity in r t 7 • Using the generalized GR ap- 
proach, we show in Appendix [C] that the expression for 
the transport scattering rate obtained in this way remains 
valid also for <5 — > 0. 



B. Transport scattering rate: results 



FIG. 12: (Color online) Transport scattering rate for a g N = 
4 x 1CP 3 (double logarithmic scale). Dots: exact values ob- 
tained by numerical evaluation; solid lines: analytical asymp- 
totics, Eq. (J5EJ). 



IV. Again, this is not so for energies that are in the regime 
II. In this case, we have an additional contribution from 
region 1, like in section llV Al The calculation outlined 
in appendix iBl yields 



T 

N 



a g T, 



a 2 g NT 



(t) 



w '(7 



i 



ii 



in 



IV 



(68) 



The characteristic momenta dominating the transport 
scattering rate as well as the relevant regions in Q-fl 
plane as shown in Table [3] The transport scattering rate 
in Regimes II and III is dominated by momenta much 
smaller than temperature, which allows us to find the 
numerical prefactors (given in Appendix B). Regimes I 
and IV is determined by momenta of order temperature, 
which does not allow us to find the numerical coefficient 
in these regimes. The analytical asymptotics are plotted 
alongside with the exact numerical result in Fig. [12J 



T tr 
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II 


III 


IV 


qVF ~ OJ 


T 


a 2 g N 2 T/e 


c 


T 


Regions 


1,2,3,4 


1,3 


3 


3 



TABLE 3: Momentum/frequency scales and regions of the 
Q-Q, plane that dominate the transport scattering rate r t 7 in 
different domains (I, II, III, and IV) of energy e. 



Similarly to the quantum scattering rate, the transport 
rate is dominated by region 3 in the regimes II, III, and 



Similarly to Sec. IIV C[ we see a strong enhancement of 
the transport scattering rate in region II. When energy 
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approaches temperature, the GR result 



T) 



a g NT 



(69) 



is reproduced. 

Comparing the transport scattering rate with the en- 
ergy relaxation rate, 



— oc In 



a g N' 



(70) 



we observe that in the limit of small a g the relaxation of 
energy due to the inelastic collisions occurs much faster 
than the velocity relaxation, r tr 3> te- The difference 
between the two rates comes from the fact that the for- 
ward scattering in graphene is strongly enhanced. The 
RPA screening suppresses the contribution of scattering 
angles smaller than a g thus regularizing the logarithmi- 
cally divergent contribution to the energy relaxation rate. 
On the other hand, the transport factor in /C tr , see Eq. 

kills the contribution of the forward scattering to 
the transport scattering rate much more efficiently than 
the RPA screening (this is the reason why the GR result 
for Tt r is parametrically correct), so that no logarithmic 
factor arises in T tr . 

The relation Eq. ([70]) justifies the hydrodynamic ap- 
proach^—: the distribution functions of electrons and 
holes equilibrate within each type of carriers much faster 
than the direction of the velocity is changed. As a re- 
sult, the distribution functions effectively depend only on 
the velocity direction. This means that the interaction- 
induced transport scattering rate entering the observ- 
ables should be averaged over the temperature window: 



de p(e) 



(r t r 1 (e)} = - 



df 1 

de Ttr(e) 



(71) 



Evaluating these integrals numerically using the exact 
expression for the RPA propagators, we find 

(T^(e))~l.l8a 2 g NT. (72) 

As we will see below, this averaged value of the transport 
scattering rate can be substituted into the Drude expres- 
sion for the conductivity, which yields the parametrically 
correct result for the collision-dominated conductivity of 
clean graphene. 



C. Collision-limited conductivity 

Above, when calculating the conductivity, we have as- 
sumed that the finite lifetime of quasiparticles is pro- 
vided by some artificially introduced broadening 5 which 



mimics disorder. The introduction of the artificial life- 
time allowed us to identify the contribution of the in- 
elastic collisions to the transport scattering rate by con- 
sidering the perturbative-in-interaction contributions to 
the conductivity. Let us now discuss the conductivity 
of clean graphene, or, more precisely, the conductivity 
of graphene in the regime when the inelastic collisions 
dominate over disorder scattering. 

The conductivity of graphene at the Dirac point is a 
rather intricate quantity. In the absence of interaction, 
transport in the ballistic limit shows remarkable pecu- 
liarities in graphene^ The interplay between vanishing 
density of states and vanishing scattering rate leads to 
the non-universal conductivity that depends on the mea- 
surement details. In particular, finite size clean graphene 
sample of the "short-and-wide" geometry shows a behav- 
ior analogous to that of a normal diffusive metal. The 
zero-T conductivity of such a setup is a = 4 x e 2 /irh. 
The same value of conductivity was predicted for an in- 
finite sample with a large but finite electron lifetime. 
A different result, a = e 2 /2h, was found in the un- 
doped graphene at a large frequency. At any non-zero 
value of the chemical potential, the conductivity of clean 
graphene is infinite. At finite temperature, energies 
within the temperature window contribute to the con- 
ductivity. This implies that the Dirac-point conductivity 
becomes infinite in the noninteracting case at any non- 
zero T. 

As has been shown in Refs. 56,57.59, the conductivity 
of clean undoped graphene becomes finite due to the in- 
elastic electron-electron collisions. The estimate for the 
collision-limited conductivity can be obtained by substi- 
tuting the typical value of interaction-induced transport 
scattering time, Eq. (|69|) . and the typical density of ther- 
mally populated states, p(T) ~ NT/vp, into the Drude 
formula, which yields 



e 2 / x 2 , . e 2 NT 
<J = — p(T)v F T tr (T) ~ -r—i 
h h vp 



1 



a 2 g NT 



e 2 1 
h a g 



Note that the explicit dependence on T, A, and vp drops 
out from this formula; however, the temperature depen- 
dence appears in Eq. (|73[) implicitly through the renor- 
malization of o 9 . A more rigorous calculation of the 
collision-limited conductivity requires the analysis of the 
kinetic equation ! 56 ' 57 Importantly, the fast energy relax- 
ation discussed above not only simplifies such an analy- 
sis, but also reduces the kinetic equation to the hydrody- 
namic modeler—. 

The consideration of Sec. IV Al which allowed us to 
find the transport scattering rate from the expression for 
the conductivity, relied on the perturbative treatment of 
the interaction. This assumes the following hierarchy of 
the energy scales: 



t^ 1 < 5 < T. 



(74) 



The first inequality implies that the broadening of the 
Green's functions is due to the artificial "disorder" rate 
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S, whereas the second inequality establishes the ballistic 
regime which allows us to neglect the dressing of inter- 
action by disorder. Since the characteristic frequency 
transfer in the transport scattering rate is of order of 
temperature, the resulting t^ 1 does not depend on S as 
long as S <C T. Furthermore, in AppendixOwe have cal- 
culated r^ 1 from the generalized GR approach for 5 = 
and reproduced the transport kernel /Ct r from Appendix 
IB] Therefore, we expect that the Drude formula, yield- 
ing Eq. (|75)) . is applicable also for r^ 1 3> 5, i.e. in the 
collision-dominated transport regime. This expectation 
is supported by the results of Refs. [56ll57ll59l 

Assuming the validity of the Drude formula for 6 <C 
r t ~ , we evaluate the conductivity using the transport 
scattering rate given by Eq. (|72|): 



Nn 



; x (e)> 



dep(e) 



df 
de 



0.58 



(75) 



This result has the same form as found in Refs. 1561)57 
within the kinetic equation approach. The only differ- 
ence is in the numerical value of the prefactor. It is worth 
noting here that the kinetic approach of Refs. 15611571 was 
based on the GR calculation of the self-energies. As we 
have seen in Sec. IV A[ the RPA result for the trans- 
port scattering rate at relevant energies e ~ T [as well as 
the energy-averaged characteristic rate, Eq. (|71)) ] has the 
same form as given by the GR. However, the proper inclu- 
sion of the finite- T changes the prefactor in the transport 
scattering rate. This explains the difference between the 
numerical prefactors between our result, Eq. (|75|). and 
the result of Refs. |56ll57t We thus see that it is impor- 
tant to include the full RPA-interaction into the kinetic- 
equation approach. This will be done elsewhere. 

Finally, let us discuss on the qualitative level the effect 
of disorder on the Dirac-point conductivity in graphene. 
For simplicity we set N = 1 below. In the absence of 
interaction, effect of disorder on transport in the ballis- 
tic regime is highly unconventional and strongly depends 
on the type of randomness i 13 ' 14 ' 22 ' 75 On the other hand, 
as we have seen above, the interaction effects are cru- 
cially important for the transport already in a clean sys- 
tem. It is thus important to explore transport in realistic 
graphene structures, with both the electron-electron in- 
teraction and disorder taken into account. 

The role of disorder here is twofold: (i) potential dis- 
order introduces velocity relaxation, thus contributing to 
the transport scattering rate: 



1 , 1 
a 2 g T + — 

Ttr Tdit 



(76) 



and (ii) establishes a finite density of states already in 
the Dirac point: 



p cx T 



1 

Tdis 



(77) 



Substituting these formulas into the Drude conductivity, 
we obtain the following result describing the crossover 




FIG. 13: Schematic plot (only parametrical scales are given) 
of the temperature dependence of the Drude conductivity (in 
units of e 2 /h; solid line) a g <C l/ln(Ardi s ), where A is the 
bandwidth. For simplicity, the logarithmic temperature cor- 
rections to a g which comes from the renormalization is not 
shown. Dashed lines: low-T behavior of the conductivity gov- 
erned by quantum interference effects; depending on the char- 
acter of disorder, the localization, antilocalization, or critical- 
ity (coincides with the solid line) may occur. Dash-dotted 
line: expected high-T behavior governed by electron-phonon 
scattering. 




T 



Ac"' 



FIG. 14: Schematic plot of the temperature dependence of the 
Drude conductivity (in units of e 2 /h) for strong interaction, 
a g ^> l/ln(Ardis)- The characteristic scale To is given by 
To ~ r^l ln 2 (A/Xb) ~ rj^ 1 ln 2 (ATdi s ). For possible deviations 
at low and high T (interference and phonon contributions, 
respectively), see Fig. 1131 



between the collision-dominated and disorder-dominated 
regimes: 



T + 



h a 2 n T+^-' 

9 Tdis 



(78) 



This expected temperature dependence of the Drude con- 
ductivity is shown in Fig. [13] for the case weak interac- 
tion (or strong disorder), a g <C l/ln(Ar ( ji s ) when the 
renormalization of a g gives only small logarithmic cor- 
rections. For stronger interaction (or weaker disorder), 
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tt g <C l/ln(Ardis)) the renormalization of 



a a 



1 



4 m T 



becomes strong [so that the renormalized coupling "for- 
gets" about its bare value, a g (T) — > 4/ln(A/T)] already 
in the collision-dominated regime, see Fig. [TU One sees 
that, since the interaction-induced transport rate con- 
tains a 2 g whereas the density of states of thermally ex- 
cited quasiparticles does not, the two crossover T-scales 
appear, which establishes an intermediate regime of the 
ballistic transport, 



1 

Tdis 



1 



Oi 2 g T Als ' 



Again, a more rigorous derivation of the conductivity 
in the presence of both disorder and inelastic scattering 
is based on kinetic-equation approach and will be per- 
formed elsewhere. Note that we neglected the phonon 
contribution to the relaxation rates (for estimate of their 
contribution, see, e.g. Ref. l24h which becomes relevant 
at sufficiently high temperatures. 

At T ~ l/idis the Drude conductivity becomes of the 
order of conductance quantum and the dephasing rate 
becomes of the order of T. At lower temperatures, the 
T dependence of the conductivity is governed by in- 
terference effects: localization, antilocalization, or crit- 
ical behavior may occur, depending on the symmetry 
of disorderJ&i^ The crossover scale l/rais in typical 
experiments on high-quality graphene is in the range 
T ~ 1 - 100-ftT. 



the inelastic scattering rates as compared to the zero- 
temperature case. We have calculated the finite-T quan- 
tum scattering rate, see Eq. (f57| and Fig. El which is 
relevant for dephasing of interference processes. We have 
identified an hierarchy of regimes, Eq. ([5]), arising due to 
the interplay of a plasmon enhancement of the scattering 
and finite-temperature screening of the interaction. The 
lifetime of quasiparticles with energies close to the Dirac 
point has been found to be independent of the coupling 
constant. We have further calculated the energy relax- 
ation rate, Eq. (|STt and Fig. [5J and transport scattering 
rate, Eq. (pS]) and Fig.JT^l For all the three rates, we have 
found a non-monotonic energy dependence which has 
been attributed to the resonant excitation of plasmons. 
Finally, we have discussed the collision-limited conduc- 
tivity of clean graphene as well as the expected behavior 
of the high-temperature conductivity in the presence of 
disorder, see Eqs. (|75|) and (175)) . respectively. Our results 
complement the kinetic-equation and hydrodynamic ap- 
proaches for the collision-limited conductivity. 

Our approach that employs the Keldysh formalism can 
be generalized for the treatment of physics of inelastic 
processes in strongly non-equilibrium setups. In partic- 
ular, this framework is expected to allow us to investi- 
gate interaction effects on full counting statistics of the 
electron transport in graphene and to develop the the- 
ory of tunneling spectroscopy in strongly biased graphene 
setup. 
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VI. CONCLUSIONS 

In conclusion, we have analyzed the inelastic electron- 
electron scattering in graphene using the Keldysh dia- 
grammatic approach. We have demonstrated that fi- 
nite temperature strongly affects the screening proper- 
ties of graphene. This, in turn, dramatically influences 
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Appendix A: Polarization Operator 



In this Appendix we evaluate the polarization operator, Eq. (|33|) . at finite temperature. In a scattering process with 
emitting a photon, three different momenta are involved which form in two dimensions triangles as shown in figure 
1151 Specifically, if an electron before scattering has the momentum p and the emitted photon carries momentum 
q, the electron that is left over has to carry p — q. The angular integration over the transferred momentum q 
becomes complicated. To proceed further it is convenient to choose elliptic coordinates defined by £ = p + |p — q| 
and r) = p — |p — q|. The corresponding coordinate system is shown in Fig. 1151 Using the elliptic coordinates, the 
expressions for the imaginary and real parts of the polarization bubble take the form 



2tt 



Re 



Q 



drj 



Q 



cosh(Q/3) + cosh(f?Q) _ p2 



cosh(Q/3) + cosh(£Q) 



(Al) 



19 



••" T| = const. 



const. 




FIG. 15: Sketch of elliptic coordinates 



RcIF* = 



Q 



sinh(£Q) 



J d V ■ 

o 



e-ip 2 -e 



sinh(77Q)i 



'£ 2 -l ^ 
1 — r\ 2 ft 2 — r\ 2 



cosh(£Q) + cosh(r7<9) 



(A2) 



Here and below we introduce the notation ft = Q/Q to decouple expansions in small or large Q from the behavior at 
the singularity Q = ft. For simplicity we set Vp = 1 and T = 1. The integrals in Eqs. (jAll) and (1A2[) are evaluated 
separately in the four regions shown in Fig. m 



a. Region 1 

In region 1 the condition Q C 1 and ft > 1 hold, which means that in Eq. (IA1[) only the part with the 77-integral 
is left. Expanding the integrands in small r/Q < 1, we get 



imn* ~ /*, 



y/l=rp _ Q tanh(Q/3/2) 



ZK^fft^lJ cosh(Q/?) + l 8 



(A3) 



for imaginary part and 



, 2 In 2 
' Q' 1 



Rell* 



-Q f 1 djq Q V 2 
* 2 Jo^fP-il 2 



d£ 



cosh(£Q) + l 



ln2/ \ft\ 



(A4) 



for the real part of H R . The term oc sinh(£Q) in Eq. (|A2I) yields a contribution of the order of ~ Q 2 . The second 
term oc sinh^Q) results in Eq. (|A4[) . 
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b. Region 2 

In region 2, we have Q ^> 1 and (3 > 1. In Eq. (|Alj) again only the 77- integral is left and by neglecting cosh^Q) 
and expanding in large Q we get 

1 



^sinMQ^) r ^ = 1 Qsig^ 

2^ Jfr^T J cosh(g/?) 8 v ; 



By expanding Eq. (IA2[) in large Q and resolving exponentials in the denominator by a geometric series, we see that the 
leading contribution to the real part of the polarization operator is small in region 2, being of the order of ~ Q~ 2 <C 1. 



c. Region 3 



Region 3 is characterized by the conditions Q < 1 and j3 < 1. In this region, only the £- integral contributes to 
Eq. (|A1[) . By expanding the integrand in small Q/3 we get 



. 2 hi 2 

' Q 2 



Iran* * A* * ^Re^. (A6) 

2^^/13^27 M+cosh(cQ) 7T <• ; 



The simplification of Eq. (| A2|) in region 3 is similar to that leading to Eq. (IA4[) . The contribution of the first term 
in Eq. (|A2|) is of order QlnQ, yielding 



, 2 In 2 



^ 2 7o ^/T^7f P 2 - V 2 J S cosh(£Q) + l 7T 



d. Region 4 

Finally, in region 4 the conditions Q > 1 and (3 < 1 are fulfilled. In Eq. (|A1[) only the ^-integral contributes and 
by neglecting cosh(Q/3) and expanding the integrand in large Q£ we get 

00 

nan* * 0= f^Ve^Te^ * ^m^L. (A8) 



1 



In region 4 the simplification of Eq. (|A2[) differs from that in region 2 in one important point. By neglecting cosh(Q?7) 
in comparison to cosh(Q£), the first term [oc sinh(Q£)] yields no principle value integral, while /3 is smaller 1. This 
results in 

00 i. 

q r„ e 



Re n« ~ df > / d^yi^ = ± . (A9) 



1 v '0 



Appendix B: Calculating the rates 



In this Appendix we calculate the integrals for inelastic scattering rates. 
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1. Definitions of the rates 



The rates we are interested in are denned by inserting the kernel JCj into the integrand for the imaginary part of 
the total self-energy, Eq. (1511) (j = q, +,E,tr): 

• Total quantum scattering-rate (j = q): 

/C q = 1; 



Energy relaxation-rate (j = E): 

r < 



Transport scattering-rate (j = tr): 



r 1 • 2 a 1 q 2 sin 2 7 
A-tr = 77 sin — 



2 2 e 2 + q 2 - 2egcos7 



Here is the angle between incoming particle and outgoing particle and 7 is the angle between q and p. The origin 
of /Ctr is explained in Sec. IV Al The kernel for the chiral scattering rate r^ 1 depends on the scattering channel. In 
particular, for + chirality at e > the electron-electron scattering kernel contains cos 2 (9/2) and the electron- hole 
scattering kernel contains sin 2 (6/2) due to Dirac factors. 

It is convenient to introduce the dimensionless energy y = e/2T. The integrals we have to handle are of the following 
form: 

7 T 1 (y) = -^/^Q/^ E [/C j (0,Q,j/)Im£>«, A (0,Q)[coth(fi)+tanh(y-0)]]. (Bl) 

F n=y±y'y 2 +Q 2 +2Qycos~f 

Since the combinations y ± \Jy 2 + Q 2 + 2Qy cos 7 lead to complicated integrands, we split the integrals into the parts 
corresponding to Q <C y and y Q. Here 7 is the angle between the transferred momentum q and the initial 
momentum p. 



2. Simplifying the integrand 



The imaginary part of the interaction propagator is given by 

Dl(Q) 

[1 + D (Q) N Re n(fi, Q)Y + [D (Q)iVIm U(Q, Q)] 



where 



Do{Q) =TQ (B3) 



is the bare Coulomb interaction. 

The splitting of the Q integral leads to the following simplification: 



2fj 



^ ±i ,-Q» 7 Q«, VW^l^l ?^T'. (B4) 
ycos7 y < Q j 2vw|sm-j| 

2VyQ I cos i I 



Here +, — correspond to the two possible values for ft that appear in Eq. (|B1[) . These signs reflect the separation 
of integration domains into the parts above and below the mass-shell line Q — fl: + corresponds to 57 > Q and — 
corresponds to < Q. Below we simplify the integrands I = [g + /] /Cj Im_Dp PA separately in each to the region in 
the Q vs £1 plane. 
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a. Region 1 

In region 1, the dominant contribution to all the rates comes from the domain y -C Q. Using the simplified 
polarization operator, Eq. (|40p . we obtain 



-ri i n \ v f 1 71-2 



Vy^ l sin i I 



T N 8 /2v?V^ 



x < 



■ 7 
sin — 
2 



In 2 



1. 



j = q 



2cos 2 ^, j = + 

4Q 2 , j = E 
l 

2 



(B5) 



1-2 

sin 7, J = tr 



Here and below the first digit in the subscript of the function X denotes the region in the Q vs plane, while the 
second digit is 1 for Q <C y and 2 for Q ^> y. 



b. Region 2 



The contribution to all the rates coming from the region 2 is exponentially small in Regime IV. In other Regimes, 
it contains at least an extra a g as compared to the contributions of region 3, except for the situations when results for 
the rate are determined by momenta q of order T/vf- In this situations, it turns out that the asymptotics produced 
by region 2 is the same as the asymptotics of region 3. We remind the reader that finding the numerical value of the 
prefactor is beyond our analytical approach when integrals are dominated by q ~ T/vf- Therefore, there is no case 
where we need to calculate the contribution of region 2. 



c. Region 3 

Region 3 appears to be the most important region because most of the final results for the rates are determined by 
this region. In this region both small and large y compared to Q are important: 

• 2/> Q : 



v < Q 



\smj\Q- 1 



NT In 2 



Q 



a g N In 2 



sin 7+1 



1, j 

2, j 
4Q 2 cos 2 7, j 
i^sin 2 7 , j 



E 

tr 



(B6) 



£3,2(2/, <3, 7) 



1 2tt 



-3/2, 



Q 



a g N In 2 



cos 



1, 




j 


= q 


2 sin 2 


2 

2 ' 


j 


= + 


4Q 2 


, 


j 


= E 


isin 2 


7, 


j 


= tr 



(B7) 



d. Region 4 

This region is only important for the energy relaxation rate (j = E). The corresponding contribution is governed 
by Q *C y, so that only the energy range y 3> 1 is of interest, where 

t-e 1 n \ VF a2 g n2/LN sin 7 cos 2 7 

£4,1 (y, Q, 7) - -?f7 /!J - n — 2 ( B8 ) 

1 ytZirQ ( a a irN\ 
' sin 7 H — ^- — 



3. Results for the rates 



Finally we estimate the rates in the Regimes I, II, III, and IV as defined in Fig. [51 using the simplified integrands 
of the previous part. 
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a. Region 1 



The contributions of Region 1 to the quantum and transport scattering rates are only relevant for energies in Regime 
II for y <C Q, where we have: 



min(y,l) 

Performing first the angular integration over 7, we obtain the quantum scattering rate: 



(B9) 



Tq 1 {y) - 4a 5 T J &a-j=L= ~ 27ra g T, 



(BIO) 



the chiral scattering rate: 



tT/ 1 (y) ~ 4a g T / da \/l — a 2 ~ Tra g T, 



(Bll) 



and the transport scattering rate: 



1 

r" 1 ^) ~ 8a g T Jdaa 2 \/l - a 2 ~ ^a g T. 



(B12) 



Here the integrals over the variable a = a g N AyQ correspond to the Q integration. 



Regime 




Q 


ft 


7 


e 


II 


~ 27ra 9 T 




<9 


< 7 < TV 


< < 7T 



TABLE 4: Contribution of region 1 to the chiral quantum scattering rate r + 1 in Regime II and the characteristic values of 
<2,fi, 7, and 6 dominating this contribution. 



Regime 




Q 




7 


e 


II 


~ na g T/2 




Q 


< 7 < TT 


< < TT 



TABLE 5: Contribution of region 1 to the transport scattering rate T tr in Regime II and the characteristic values of 7, 
and 9 dominating this contribution. 



In all the regimes, the contribution of region 1 to the energy relaxation rate is parametrically the same as that of 
other regions. Since the numerical prefactor is not accessible within our calculation (the integrals are dominated by 
qvF ~ T), we have chosen to present only the calculation of the contribution of region 3. 



b. Region 2 

There is no important contribution from region 2 (see Sec. IB 2 bl above). 
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c. Region 3 



This is the most important region for all the rates in most of the Regimes. The corresponding integrals are expressed 
through the kernels I 3 3 4 introduced in Sec. IB 21 as follows: 



miiXy- 1 ) 27T 




2T 2 



Using the short-hand notations 



1 



a g N In 2 ' 



a g Nlii2' 



we obtain the following result for the contribution of region 3 to the quantum scattering rate 

linfx' ,y' 



-T< 



arsinh 



dx- 



y/x{x + 2) 



^/x(x + 2)(x + 1) 



arsinh 



dx- 



s/') (i + ^)yi + 4y'(v^+^ 



(B13) 



(B14) 



Regime 




Q 


Q 


7 


9 


I 


jv J 


1 


-Q 


< 7 < 7T 


< < TT 


II 






-Q 


< 7 < ty 


< < TT 


III 




OgNT 


Q cos 7 


< 7 < TT 




IV 


7r y 


a g NT 


Q cos 7 


< 7 < tt 





TABLE 6: Contribution of region 3 to the total quantum scattering rate t^ 1 in all energy domains and the corresponding 
characteristic values of the transferred momentum Q, transferred frequency fi, the angle 7 between momenta p and q, and the 
scattering angle 9. 



For the contribution of region 3 to the chiral scattering rate we get: 



-T< 



arsinh 



dx- 



^x{x + 2) 



sjx{x + 2)(x + 1) 



da; 

minfx' ,y') 



1 + 4y'(Vx~ + ^) arsinh [v^V^ + rj=)] - VW(Vx + 7;) 



(1 + x) 



2 4y' 



(B15) 



Regime 




Q 




7 


9 


I 


jv 


1 


-Q 


< 7 < 7T 


< <9 < TT 


II 


~ TTOgT 


oigN 2 /y 


-Q 


< 7 < TT 


<9 <TT 


III 




a g NT 


Q cos 7 


< 7 < TT 


< ^ « 1 


IV 


7r y 


a g NT 


Q cos 7 


< 7 < TT 





TABLE 7: Contribution of region 3 to the chiral quantum scattering rate r, 1 in all energy domains and the corresponding 
characteristic values of Q, £1, 7, and 9. 
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The contribution of region 3 to the energy relaxation rate reads: 

iin(V,i/) 

,vHv)*-%*tI f dx- 



ir{x') 2 



arsinh ( \/ x(x + 2)^ 
(x + I)" l yjx{x + 2) 



x 1 arsinh 



dx- 



i 



mta(x',y') (1 + x) W 1 + Ay' (y/x + 



(B16) 



Regime 




Q 




7 


6» 


I 


JV J 


l 


-Q 


< 7 < ty 


< 6» < 7T 


II 


„2rp JV i 




) 


l 


-Q 


TV ^ 






III 


9 sjy 




) 


l 


-Q 


7T ^=r- 






IV 


~ o-g AT In | 






l 


-Q 


7T ! W- 


0< 2a* < i 





TABLE 8: Contribution of region 3 to the energy relaxation rate t^ 1 obtained in all energy domains within the energy-diffusion 
approximation and the corresponding characteristic values of Q, £1, 7, and 9. Since in all energy domains the result is dominated 
by Q ~ I , the numerical prefactors can not be obtained from the asymptotics at Q <C 1 and Q ^> 1. 



Finally, region 3 yields the following result for r tr 1 : 



2ixy 



11 - 



T dx- 



x + 2 



1 - 



arsinh I yj x{x + 2) 



(a; + 1)^1(0: + 2) 



—a T dx\ — 



min(x' ,y') 



z 2 + 3 
3z 4 



arsinhz 



z 5 (l + z 2 ) 



-1/2 



(B17) 



Regime 


r tr 1 {y) 


Q 




7 


9 


I 


~ — T 


1 


-0 


< 7 < 7T 


< 9 < TV 


II 


- f^T 


a 


-Q 


CgN 




III 


y 9 


2/ 


2/(1 -2| sin 2|) 


v v 


ff<^«l 

y 


IV 




1 


QC0S7 


y 


6/ < fVY. < 1 



TABLE 9: Contribution of region 3 to the transport scattering rate r^ 1 in all energy domains and the corresponding charac- 
teristic values of Q, Q, 7, and 9. Since in energy domains I and IV the result is dominated by Q ~ 1, the numerical prefactors 
in these domains are beyond the accuracy of our approximations. Furthermore, since in Regime III the integral is dominated 
by Q ~ y, the splitting of the integrand according Eqs. (|B6|) and (|B7|) does not reproduce the correct prefactor. 



d. Region 4 



This region is only relevant for the energy relaxation rate in Regime IV, where the main contribution comes from 
Q < V, yielding 

max(i/,l) f - "5T 

_ x 8T 2 f dQ f d7 T E, 7 «f-i 127V 2 3/2 / 16 \ 

Te(2/) = ^T J 2t7 Q J 2^ ff (»,0,7)«^o ir r» /ln^__j. (B18) 

1 

Since the result is determined by Q ~ ~ y 3> 1, the energy-diffusion model is, in fact, not applicable in this Regime. 
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FIG. 16: Diagram for the "elementary" inelastic scattering amplitude. 



Appendix C: Scattering Rates from the Generalized Golden Rule 

In this Appendix we analyze the golden rule approach to calculating quantum and transport scattering rates. The 
quantum scattering rate for electrons was introduced in Eq. (|52l) . Using the imaginary part of the self energy Eq. 
(24), we obtain 

T+ 1 (e)oc J dq J du; \D KPA (q, w)| 2 [1 - n F (e - w)] J dpTr [V + (p) S(e - v F p) Ao(e - u, p - q)] 

x J dE [n F (E - lj/2) - n F (E + lj/2)] J dkTr [A (E - lj/2, k - q/2)A {E + lj/2, k + q/2)] . (CI) 

The spectral weights Aq are defined in Eq. ([5]). Equation (jCll) is completely equivalent to Eq. (46) of the main text. 

We can regard the quantum scattering rate as the probability (per unit time) of the electron decay with emitting 
an electron- hole pair. The amplitude of this decay process is given by the "half of the self-energy diagram shown 
in Fig. [THJ The incoming and outgoing particles in this diagram are taken at the mass shell. Equation (|C1|) has the 
form of the Fermi golden rule with the amplitude determined by the RPA-screened interaction. 

We can now apply the golden rule to calculate the transport scattering rate. This amounts to including an extra 
transport factor [v, • (v; — Vf)]/v F accounting for the change of the current due to scattering, in the integrand of Eq. 
(jCip . Here Vj and v/ are total velocities of incoming and outgoing particles. For the linear electronic dispersion the 
velocities of individual quasiparticles are determined by the relation 



de _ vpp 
dp e 



(C2) 



This results in the following transport factor: 



Vi-Cvi-Vf) _ iV(P P q k+q/2 , k-q/2 \ 

v* - Vf 7{7 ~eT^]2 + e - lo/2 ) ' 

Note that in the conventional case of massive particles with the quadratic electronic dispersion, v = p/m and hence 
Vj — v/ cx p — (p — q) — (k + q/2) + (k — q/2) = 0, which reflects the fact that because of the total momentum 
conservation there is no current relaxation in conventional metals due to the electron-electron interaction. 

We now apply the particle-hole symmetry of the graphene spectrum in order to simplify the expression for transport 
scattering rate. Let us reverse the integration variables in the second line of Eq. (jCll) : (E, k) n- (— E, — k). Using 
the symmetry A^e^p) = Ao{— e, — p), we see that the trace in the integrand is not changed. The difference of 
two equilibrium distribution functions in the second line of Eq. (|C1[) is also independent of the sign of E. Thus 
the contribution of particle-hole pair into the current relaxation vanishes, which corresponds to the absence of the 
Coulomb-drag contribution to the conductivity^ 7 in undoped graphene (Fig. II Of ). This allows us to keep only the 
part of the transport factor Eq. (|C3|) which is even under reversing E and k, 



2 p / p p- q 

f P - 



= 1 - sign(e - oj) cos 9, (C4) 
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where 9 is the scattering angle. Inserting this reduced transport factor into the first line of Eq. (|C1|) and calculating 
the trace of projection operators, we obtain 

[1 - sign(e - w) cos 9} Tr [V+ (p)Ao(e u, p q)] 

= — [l — sign(e — uj) cosfl] [(1 + cos8)S(e — uj — vp\p — q\) + (1 — cos8)S(e — ui + vp\p — q|)] 

= i(l-cos 2 0)[(5( e - W -^|p-q|)+ ( 5(e- W + ^|p-q|)]. (C5) 

This way the transport factor 1 — cos 2 9 appears for Dirac fermions. 

The golden rule calculation of the transport scattering rate reproduces the result obtained from the Drude conduc- 
tivity Eq. (|64|) . Note that an additional broadening 6 was not used in this calculation. This shows that our result for 
can be actually applied in the limit t^ 1 ^> S as we argued in Sec. IV CI 
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